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Abstract 

Recently the issue of EPR-like correlations in the mutual probability of detecting 
neutrino together with accompanying charged lepton has received a new impetus. In 
this paper we describe this effect using the propagators of the particles involved in 
the Schwinger’s parametric integral representation. We find this description more 
simple and more suitable to the purpose than the usual momentum-space analysis. We 
consider the cases of monochromatic neutrino source, wave packet source, and neutrino 
creation in a localised space-time region. In the latter case we note that the space-time 
oscillation amplitude depends on the values of the neutrino masses, and becomes rather 
small for large relative mass differences (mass hierarchy). We obtain the expressions 
for the oscillation and coherence lengths in various circumstances. In the region of 
overlap our results confirm those of Dolgov et al. Q 
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1 Introduction 


Space-time oscillation of neutrino flavour is considered to be the most promising effect 
observation of which might indirectly establish nonzero neutrino mass. By its very nature it 
requires spatio-temporal description of the processes of neutrino creation, propagation and 
detection, and of the similar processes that occur with the accompanying particles. Such 
a description has been performed in and further developed in |^, Q without ambiguities 
that sometimes accompany noncritical use of neutrino flavour eigenstates. 

Recently BB the issue of EPR-like correlations in the mutual probability of detecting 
both neutrino and an accompanying charged lepton has received a new impetus. (It was 
previously considered in |^.) In order to simplify derivation of the basic effects the authors 
of m combined together descriptions in the configuration space and in the momentum space, 
of the same relevant processes, using simultaneously such mutually exclusive notions as 
sharp wave packets in momentum space and dehnite space-time a posteriori trajectories 
of particles. EPR-like experiments of the same type involving neutral kaon and B meson 
oscillations were considered in p. In this paper the authors also adopted a simple approach 
using the action values on the particle classical trajectories to evaluate the relevant phase 
factors in the probability amplitude. Although the results obtained in such a simplified 
approach are correct, they also might call for a more careful derivation. This will be the aim 
of the present paper. 

In this paper we try to analyse the phenomenon in a consistent way using the propagators 
of the particles involved in the Schwinger’s parametric integral representation (see Eq. (0) 
below). We hnd this description more simple and more suitable to the purpose than the 
usual one which employs propagators in momentum space representation. This latter in¬ 
volves rather complicated momentum integrations (see, e.g., p, and, it seems, frequently 
obscures the physical picture of the phenomenon. Our treatment will be general and will 
contain the analysis of the EPR-like experiments of detecting neutrino together with the 
accompanying charged lepton, as well as the standard textbook examples of neutrino flavour 
space-time oscillations. 

After preliminaries in the following section, in Sec. ^ we consider the case of a monochro¬ 
matic neutrino source and the probability of mutual detection of the neutrino and of the 
accompanying charged lepton. In Sec. ^the effect of wave packet neutrino source is analysed. 
We obtain the expressions for the oscillation and coherence lengths in various circumstances. 
The case of neutrino source in a strongly localised space-time region will then be considered 
in Sec. In this case the space-time oscillation amplitude depends rather strongly on the 
values of the neutrino masses, and becomes rather small for large relative mass differences 
(neutrino mass hierarchy). We summarise our results in Sec. ^ In the Appendix we provide 
an alternative derivation of the probability amplitude for the case of monochromatic neu¬ 
trino source, in order to elucidate the difference between this case and the case of neutrino 
source strongly localised in space-time. 
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2 Preliminaries 


Throughout this paper we consider a process in which neutrino is created together with 
accompanying charged lepton, and afterwards both particles are detected. The charged weak 
currents lOai^ are involved in the description of this process, where Oa = 7q-( 1 + 75), and 
75 = *7 o7i 7273- The amplitude of creation of a l-u pair at space-time point x is proportional 
to OaJsi^), with Js{x) being the source current responsible for this process. In the case of 
pion we would have Js{x) oc d°‘(j)n{x), where (pnix) is the pion wave function. The charged 
lepton produced at space-time point Xc in a flavour state a can propagate to space-time point 
xi, and the neutrino to space-time point x„, at which points these particles may be detected. 
At the space-time point Xn one may detect neutrino-induced charged lepton production 
of flavour b. The amplitude of such a process in which neutrino and the corresponding 
antilepton are created and subsequently detected will contain the factor 

OpJnY^UljUja j dXcSj{Xn,Xc)OaJs{Xc)Sa{Xc,Xl) , (1) 

j 

where Uja is the unitary matrix of neutrino mass-flavour mixing amplitudes, Ulj is its Her- 
mitean conjugate, a and b numerate flavours, j numerates the neutrino mass eigenstates, Sj 
and Sa are, correspondingly, the Feynman propagators of neutrino mass specie j with mass 
rrij and of charged lepton flavour a, and is the current involved in the neutrino detection 
process, localised around the space-time point Xn- 

The total amplitude that will describe the detection of charged lepton and neutrino 
events will contain, besides the factor (0), also positive and negative energy wave functions 
of different hnite particles involved in the detection process. These factors are of particular 
nature, they do not affect the dependence of the amplitude on the space-time coordinates 
Xi and Xn, hence they will be omitted as irrelevant to the main topic of this paper. Due to 
these factors, however, as well as, in typical cases, due to the positive-frequency character 
of the source, the integration region over Xc in (|^ will be effectively restricted to the causal 
past of both points Xn and x;. 

The Feynman propagator S{x, y) = Sm{x — y) for the Dirac held of mass m has the form 

Sm{x) = (i7“^a + m) Dm{x ), (2) 


where Dm{x) is the Feynman propagator for the Klein-Gordon held of mass m. This propa¬ 
gator has the parametric integral representation (hrst considered by Schwinger, Dyson and 


Feynman in the papers collected in |]I 0 | ) 


Dm{x') 


-r lim 

Svr^ e^+O. 


dAexp 


i 

^\x^ + \ 

— ie 

)1 

~2 

\ A 


)\ 


( 3 ) 


where x^ = x ■ x = x“Xq, is the Lorentz interval squared. The factors of type in the 
amplitude Eq. (|ID will have an ehect that in the neutrino propagator the term proportional 
to a unit matrix will not contribute, and only that proportional to the Dirac gamma matrices 
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will remain. This general property is due to the equality OaOp = 0 . Thus in (J^) we can 
replace Sj{xn,Xc) by Sj{xn,x^) = Sm^ixn - xj, where 


Sm{,x) = . 


( 4 ) 


3 Monochromatic source 

In this section we investigate the case of a monochromatic source current Js{x) that can 
arise, for instance, in the process of pion decay. Let 


J^{x) oc 


( 5 ) 


with constant four-momentum p. In the case of pion we would have J^ix) oc oc 

p°‘ exp{—ip ■ x), where (pnix) is the pion wave function. We make notation 


( 6 ) 



In the amplitude (|ID we represent the propagators using (@)-(@), first perform integration over 
Xc, then over the parameters Xi and Xn that appear in the representation (|]), respectively, 
for charged lepton and neutrino propagators. The integral over Xc is Gaussian, hence it can 
be evaluated exactly, pre-exponential factors can be obtained after integration over Xc by 
taking partial derivatives with respect to xi and Xn according to (^. The remaining integral 
over A/, A„ will be afterwards evaluated in the stationary phase approximation. 

Consider the integral over Xc of one of the terms in the sum of (P. The phase in the 
exponent of the integrand will stem from the expression (^) for propagators, and from the 
source current in (|I]). It will be given by 


0 = - - 2 ^nxl^ -p-Xc. 


( 7 ) 


Its extremal point Xc = Xc{Xi, Xn) is determined from the equation 



( 8 ) 


We also have for the matrix of the second derivatives 



( 9 ) 


so that integration over Xc will produce a factor 

J . o 



( 10 ) 


where 0* is the value of the phase 0 at the extremal point: 


XlXnXh + 2p- {XlXl + XnXn) 

2 (A; -|- A„) 


( 11 ) 
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Now consider the integral over the A’s. It will be evalnated in the stationary phase 
approximation. The phase of the integrand is given by 


$ 


m. 


mt 


A, ^ A. 


( 12 ) 


where mi and are the masses, respectively, of charged lepton and of nentrino. The sta¬ 
tionary point is determined by differentiating (|T^ nsing (pT]), or by the eqnivalent conditions 
in the convenient form obtained nsing and (H), 


d\i~ 2 2A2 ’ dXn~ 2 ^^ 2X1 


(13) 


In these eqnations Xc = Xc{Xi, A„) is the solntion of Eq. (||). From (^) we have the relation 


XlXlc Pi , XnXnc Pn ; 


(14) 


satished by the extremal valnes of A’s, where pi and are the fonr-momenta that the 
charged lepton and the nentrino respectively wonld have were they free classical particles 
moving from the space-time creation point Xc respectively to the registration points xi and 
Xn- Then Eq. (^ expresses the energy-momentnm conservation law, the condition from 
which the extremal point Xc with extremal A’s can be fonnd most easily. 


We need also the matrix of the second derivatives of <h over A’s at the extremal point. 
Differentiating the identity (H) we hnd 


dXc _ Xlc dXc _ Xnc 

dXi Xi + Xn' dXn Xi + Xn' 


and, differentiating (|^, 


^2$ ^ 

W 

^ 2 $ ^ 

^ 2 $ 

dXidXn 


2 2 
Xfc _ m 

Xl + Xn Xf 


Xi + A„ A^ 

^Ic * ^nc _ 

Xl -|- Xn 


-mf 

-ml 


2Xl Xn 

+ ^n) 

2Xn + Xl 

+ -^n-) 


Pi -Pn 

XlXn (A; -|- A^) 


(15) 


(16) 

(17) 

(18) 


Let 

m = y/p ■ p , (19) 

be the effective mass of the sonrce. In the case of pion decay this will be eqnal to the pion 
mass TTijr. In a realistic case 


m — m; S> m„, mi 3> . 


( 20 ) 


Below we will see [cf. Eq. (^)] that in the limit > 0 the extremal valnes of A’s remain 
hnite. Thns we can approximate the determinant of the matrix /dXidXj by its limit as 
mn —>■ 0. The resnlt is 


/ A / y 


Pi ■ Pn 


XlXn {Xl -|- A„) 


( 21 ) 
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Therefore the integral over A’s will produce a factor 


(A/ T A 72 ) 


( 22 ) 


PrPn 

where is the extremal value of the phase <h, which is given by 

rnn\[^c-P ■ (23) 

with Xc being the extremal point of 0 at extremal values of A’s - the solution to (|[), (|T^) . 
Using the extremality conditions (||), (|T3|) we easily find 


dxi 


= -Pu 


dXr, 


= -Pr 


and 


= -Pi- Xi-Pn- Xn= -P-Xi-Pn-Xnl- 


(24) 


(25) 


Note that the four-momenta pi and pn lie in the plane formed by the four-vectors p and Xnu 
and are determined by energy-momentum conservation. 

Combining together the factors calculated m and (^), dropping the resulting overall 
numerical constant i/Sn, and using (|25|) we obtain the expression for the amplitude (|l|) in 
the case of monochromatic source 


-ip-xi ^ 


AzA 


I'^n 


(Az -h Xn)pi -P: 


- OolaPnOs (mi - UljUjaC 


—ipn-Xnl 


(26) 


where Od = OaJ^^ Os = OaJs- Note that the coordinate-dependence of the current Js{,x) 
has transformed to the phase of (^). Also note that the extremal values of A’s as well as the 
four-momenta pi and pn under the sum (^Bf ) depend on the neutrino specie j. However, the 
prefactors in our expression (^), as well as in (|30D below, are calculated only up to terms 
proportional to m^, with this precision they can be taken in the limit = 0 . 


The extremal values of Az and A„ determined by the system of equations (|]), (|^) can 
be easily obtained from the kinematics of the problem. Let us denote by t and by d corre¬ 
spondingly the time difference and the absolute spatial distance between the events Xn and 
Xi in the rest frame of the source [in which p" = (m, 0)], and by vi and the velocities, 
respectively, of the charged lepton and of the neutrino in this frame (see Fig. |^. In the rest 
frame of the source one has 


tic 


d-Vnt , I 

) di . \xiq 

Vl+Vn 


Xltlc ) 


d -|- vit 

Vl+Vn' 


dn ■ 


Xr 


Xfitnc ) 


(27) 


where also di denotes the spatial distance in the source rest frame between the point xi 
of charged lepton detection and the extremal point Xc, and dn has the same meaning for 
neutrino. Then 



( 28 ) 
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and, using ([T3|), we obtain 


Az — 



viEi 

di 



(29) 


where Ei and En are the energies, respectively, of the charged lepton and of the neutrino in 
the rest frame of the source. In this notation and in the approximation of = 0 for the 
prefactors (but not for the phase) the amplitude (PBf) will acquire the form 


g-ip-xi 

md 


Oo'^aPnOs 


3 


(30) 


By the way, from the expressions (P9|) it is clear that the extremal values of A’s remain hnite 
in the limit of rrin ^ 0, as was stated above. 



Figure 1: Time and length dehnitions in the neutrino source rest frame. 


We shall now estimate the applicability limits of the stationary phase approximation 
used. Our approximation will be good when the extremal values of A’s are much larger than 
their dispersions determined by the matrix /dXidXj. Using (p9|) and ([T^)-([T8D we obtain 
after straightforward analysis the conditions 

di < dn mdf or dn < di md^ , (31) 

under which our approximation is valid. They imply also the condition 

md S> 1, (32) 


which is quite reasonable. 

In the limit of (^) the energy-momenta pi and change relatively very slightly with the 
neutrino mass specie j, and, we remember, our prefactors in ( |30[) were actually calculated 
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in this limit. In this case the space-time behaviour of the l-u pair detection probability is 
given by 


PbaiXn,Xl) OC — 






jaUak^kb COS (^Pjk ' Xnl + 'pjk) 


(33) 


where pj and pk are neutrino four-momenta of mass species, correspondingly, j and k, pjk = 
Pj —pk, and are constant phases that stem from the product of matrices U and t/b The 
second sum in the square brackets of (|3^) describes space-time oscillations of the probability. 


From Eq. one obtains the oscillation length and oscillation time of the probability 
considered. Using the energy-momentum conservation in the source rest frame of reference 
one has 


0 


\p\jk = \Pj\ - \Pk\ 


^jk 




2vim 






2Vr,mEr, 


(34) 


where approximation uses the assumption (|20|) that neutrino has very small mass, and 


- ml. 


(35) 


Thus, oscillation length Lose and oscillation time T^sc of the (j7c)-component of (^) in this 
frame of reference are given, respectively, by (we use the limit of Vn = 1) 


_ — IUlt 

-^OSC 771 OSC jp ^ ! 


(36) 


where 

T - 

J-k — T-T- 

is the standard expression. To proceed to any other reference frame what one has to do 
is to transform the four-vector components p^j^ obtained in (|3^) to this new frame. The 
expressions (PB|) and the relevant expressions in the laboratory frame of reference have been 
obtained in [jl[. 



If one of the particles, charged lepton or neutrino, is not observed, then the probability 
of detecting the other one is uniform in space and time. This is quite obvious without any 
calculations and is due to the fact that the l-v pair creation probability for a monochromatic 
source is homogeneous in space and time. If neutrino is not detected, oscillations in the 
charged lepton detection probability disappear also due to orthogonality of neutrino mass 
eigenstates. This last cause will operate with any type of source, not necessarily monochro¬ 
matic. Specifically, it is the necessity of summing the probability over the neutrino flavour 
index b that will eliminate the oscillatory terms in this case. A detailed discussion of these 
issues is presented in [Q. 


4 Wave-packet source 

First of all let us analyse in a little more detail the effective region of integration over Xc in 
(2) in the case of monochromatic source considered in the previous section. In other words, it 
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is the region of constructive interference, from which most of the contribution to the integral 
in (|l|) comes. The extension of this region in space-time around the extremal point Xc is 
determined by the covariance matrix (|^) for given values of A’s, and by the variation of A’s 
that are determined by the covariance matrix 8“^^/dXidXj with the components (|^-(^). 
First, using (|^) and (pQ]) we obtain the estimate of the linear dimensions 5x of the effective 
region of integration for hxed extremal values of A’s as 


&X (A; + A„) 



(38) 


Next, we must estimate the linear dimensions 5xc of the spread of the extremal value 
Xc{Xi, Xn) caused by the spread 6X of the values of A’s. This latter spread can be estimated 
using ( 0 ) as 

- 1/4 


hA ~ 


det 


'_d^' 

dXidXj ^ 


\ 


XlXn{Xl -|- A,; 

PrPn 


Then using (pISD , 
the estimate 


(1^) and the condition hA A;^„ provided by Eq. 
dSX 


" A, + A, 


~ d. 


A/A 


I'^n 


Pi -Pni^l + A. 


(39) 

we will have 

(40) 


and d^) use smallness of the neutrino mass and become equal- 


Approximations in ([ 
ities in the limit —a 0. The last value in (|38|) is larger than that in (^OD hence the 
dimension of the region of constructive interference will be estimated by 


In realistic situations the source of neutrinos can often be approximated by a wave packet 
with sharp energy and momentum distribution (therefore small relative energy-momentum 
spread). Let us suppose that in the rest frame of the source the spread in the coordinate 
space is in momentum space dp, so that cr^cTp ~ 1. The source has also hnite coherence 
time, in the case of pion this will be determined by its lifetime, or by its collision time with 
the environment. Typically, however, the time spread at of the wave packet is much larger 
than its spatial spread and the latter will determine most of the interesting effects. 


If the source is to a high precision monochromatic so that its spatial and temporal spread 
is sufficiently large, namely, if 


6x -C da;, at , 


(41) 


then we can use the expressions from the previous section for the detection probability ampli¬ 
tude whenever the effective region of integration over Xc (region of constructive interference) 
lies well within the source wave packet. Due to Eq. ( |38[) the conditions (^T]) essentially imply 


d ^ En X min 




(42) 


Spatio-temporal oscillations in the mutual detection probability can be observed only up 
to certain relative distances between the detection points of charged lepton and neutrino. 
We are going to determine such maximal distances, called coherence lengths beyond which 
oscillations cease to occur. The reason for such distances to exist is that for different neutrino 
masses ruj and ruk the corresponding centres (extremal points) Xj and Xk of the integration 
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region in Xc are different. If they become sufficiently separated in space-time they may no 
longer be able to lie simultaneously within the wave packet of the source, thus components 
in the probability amplitude that correspond to different neutrino mass species will not be 
able to interfere. 

Consider this effect quantitatively. The shift four-vector Xjk = Xj — Xk lies in the plane 
of the four-momenta pi and and can be decomposed into components Xjl and Xj^ that 
go, respectively, along the directions of pi and These components can be easily esti¬ 
mated. Using Eq. (^7]) we obtain for the time components (approximation uses ~ 1 and 
Avu An„ -C vu Vn) 


fii) , 

^jk 

An) 


AVr, 


{vi + Vr, 


{d -1- Vlt) 


AjkEi 

2mE‘i 


dfi 1 




{vi + V, 


2m?El 


where Ajk is given by (^). The distances di and d^ then change as 


Adi = -Adn = Vn t^^k - '^i i'll ~ 


Ajkmf 
^ 2m?El 


Ai) 


di -\- 


A, 


jk 


2mEr, 


dn 


(43) 

(44) 


(45) 


For large enough absolute values of tjk = +4^^ and/or Adi the centres Xj and Xk will not 

be able to he simultaneously within the source wave packet. Components in the probability 
amplitude which correspond to different neutrino mass species will not be able to interfere, 
and probability will cease to oscillate. This will determine coherence lengths of the detection 
probability oscillations in various cases. 


In the case 


dn. 


mE„ 


■ di 


(46) 


from (^51 ) it follows that Adi = 0, and the shift Xjk is in the temporal direction in the source 
rest frame. For sufficiently large value oi d = di + dn the absolute value of the shift tjk 
becomes larger than the extension at of the source wave packet in the temporal direction. 
This determines coherence length - the largest value of d at which oscillations can be 
observed. In the case considered (|46D it is determined by using (|43|) , (|4^) and (^6]) as 


L, = atEnL (^1 + , (47) 

with L given by (|37|) . The condition (|4^) necessary for our approximation, will imply that 
the formula ( 0 ) is valid for 

at,— :>L 

o-t 


1 + 


mEn 


mi 


(48) 


Equation (M) implies rather special experimental coincidence conditions. In a less special 


case 


dn A 


mEn 


(49) 
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the second term in the last expression of (^51) dominates. In this case the coherence length 
will be determined by the condition that either \tjk\ becomes larger than at, or/and \Adi\ 
becomes larger than ax- Using (^31) and (|^) we obtain in this case 


T T -I 

Lc = Lm X mm [ax, at — 

El 


(50) 


The condition (|^ will determine the validity limit of (pOD to be [see Eq. 


(Jx 




> 


m 

En 


L = T^^ 


a. 


m 


(Tt, — :>—L = h 


Ct 


El 


(51) 


The case 


dn ^ 


mr 


mE„ 


■ di 


(62) 


can describe the situation when the momentum of the charged lepton is measured to a good 
accuracy by measuring its time of flight. Its analysis is quite similar to that of the previous 
cases. The oscillations in the probability will disappear at the length 


Lc = 


m?Er, 


mf 


L X mm{ax, (Jt), 


(53) 


and the analysis is valid as long as the estimate (^) is fulhlled, what gives 


(Jt (Jx mf 


(54) 


If the source wave packet is sufficiently broad in space and if one of the particles, charged 
lepton or neutrino, is not observed, the probability of observing the other one will not 
oscillate in space and time. This is because after integrating the probability (p3D over one of 
the variables {xi, Xn} the oscillatory terms are averaged to approximately zero. However, 
with the source wave packet sufficiently narrow in space (but still such that the condition 
(PI) or its equivalent (|^) holds), namely. 


(To 


•C L, 


( 55 ) 


even if the accompanying charged lepton is not observed the neutrino flavour oscillations can 
be observed relative to the source spatial position. Indeed, when integrating the detection 
probability Pi,a{xn,xi) over Xi we notice that only in a small region of Xi the probability 
will be nonzero and will be given by (|33|). This will be the region for which the extremal 
values of Xc lie within the wave packet of the source. The linear dimensions of this region are 
similar to axm/Ei, and Eq. (^) follows from the condition that these dimensions be much 
smaller than the oscillation length Lose which is given by ([36|). Then, after integration over 


Xi the phases of the oscillatory terms in (|33|) will be hxed and their dependence on the value 
of Xn will remain. Detection of neutrino alone is the case most frequently discussed in the 
literature. 
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Consider this effect more thoroughly. The phase pjk ■ Xni in the oscillating term in the 
probability (|33|) can be written in terms of the distances di and dn introduced in (|^) . 
Choosing the z-a.xis in the direction of in the source rest frame one will have 


Pjk = (p°t, 0,0, IpIjj) , = ((, 0,0, d) . 


(56) 


Taking into account the values = Ajk/2m, \p\jk ~ —^jkEi/2vnTnEn [see Eq. (|3^)] and 
using one obtains the standard expression for the phase 


Pjk ‘ ^nl 


2vnm 


dn — dl 


^jkEl 

2VnmEn 


{dn + dl) 




(57) 


From this expression it is again clear that if the distance dn is fixed with accuracy better 
than L by the position of the source wave packet relative to the neutrino detector, the 
phases of the probability oscillations will remain fixed even after integration of (|53| ) over the 
unobserved point xi, and the oscillations will be observed with respect to the value of dn- 
This condition again leads to the estimate (p5|). 

Remarkably, the phase (|57| ) does not depend on di- This fact can also be explained as 
follows. If dn is fixed and di is changing, this means that the four-vector Xni changes (say, by 
amount Axni) in the direction of the charged lepton’s four-momentum pi (this is clear from 
Fig.0). Then the oscillation phase change is pjk- Axni oc Pjk-Pi = -Api-pi = — A {pf) /2 = 0. 
The expression (|^) for the phase coincides with that derived in [Q] (see Eqs. (47) and (54) 
of this reference). 


As noted already at the end of the previous section, and as it was discussed in [|l], if 
neutrino is not detected oscillations in the charged lepton detection probability disappear in 
any case due to orthogonality of neutrino mass eigenstates. Specifically, it is the summation 
of (^^ over the neutrino flavour index b that will eliminate the oscillatory terms. 


5 Source in a localised space-time region 

First consider a hypothetical process in which neutrino, together with a charged lepton, is 
created at a fixed space-time point Xc- Note that in this case the energy and momentum 
of the neutrino created is totally undetermined. The probability amplitude of detecting 
neutrino-induced charged lepton production of flavour b at the space-time point Xn will 
contain the factor 

-^baiXmXc) j [ UhjUjaSj {Xnc) ; (^^) 

j 

if the charged lepton created together with the neutrino at point Xc is of flavour a. 

Again, as was already discussed in Sec. |^, the total amplitude that will describe the detec¬ 
tion of neutrino event will contain, besides the factor (|58D , also those related to the processes 
of creation, propagation and detection of other particles involved. These factors, however, 
are of particular nature, and do not affect the dependence of the probability amplitude of 
neutrino detection on the space-time points and Xc, hence they will be omitted. 
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The propagator S'm(x) of Eq. (|^) has the leading asymptotic behaviour (see. e.g., H) 

/ g37ri/4 \ 3/2 / ^ / /_\ /_ 

(-' ™V^) . for » 1. (59) 

~ < 1 , (60) 

271^ [x^) 

Hence, oscillations in the neutrino detection probabilities can develop in space and time 
around Xn only when rrij > 1 at least for the largest of the neutrino masses, since 

Sj{xnc) do not differ for different j in the opposite limit rrij ^x‘^^ -C 1. Consider, therefore, 

the case of 1 for all j. In this limit the amplitude (|55|) up to one and the same 

factor will be given by 

Ai,a{Xn,Xc) CX ^ Gxp . (61) 


Note the mass-dependence of the coefficients in the last equation. The space-time variation 
of the probabilities of the corresponding processes will be given by 



where 

rrijk = rrij - rrik , (63) 

and (pjl are constant phases that stem from the product of matrices U and f/h The last 
term in the equation ( |6^ describes space-time oscillations of the probabilities. 

As an illustration consider mixing between two mass eigenstates ui and U 2 , with two 
flavour eigenstates and Ue- Let, for dehniteness, mi > m 2 . 


= z^i cos 6 + 1^2 sin 9 , Ve = sin 6 + V 2 cos 9 , 


(64) 


and let neutrino be created in a flavour eigenstate The corresponding probabilities of 
detecting muon and electron events at Xn will be 


F^(x„, Xc) (X ml cos^ 9 + ml sin^ 9 + 2 (mim 2 )^^^ sin^ 9 cos^ 9 cos [ Amy 


Pe{Xn, Xc) OC 


ml + ml — 2 (mim 2 )^'^^ cos 


sin^ 9 cos^ 9 , 


(65) 

( 66 ) 


where Am = mi — m 2 . The last terms in the square brackets in the equations (|65D and 
(|66|) describe space-time oscillations of the probabilities. Due to mass-dependence of the 
coefficients in these expressions the amplitude of oscillations will be suppressed if m 2 -C mi, 
and cos6' is not too small. 


The conditions mj x^^ S> 1, together with the asymptotic form (|59|) of the propagator, 
imply that all the mass eigenstate neutrinos uj arrive at the space-time point Xn practically 
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as on-mass-shell particles, with four-velocity u = Xnd\Jdfic- opposite limit my/x^ <C 1 

independent of mass. Then, for 

2 


the Feynman propagator has the asymptotic behaviour 

instance, in our example of two neutrinos, in the region m 2 \Jxd <C 1 -C mi 
exists in the case of large relative mass difference, m 2 -C mi, one obtains 


which 


Xc) oc sin"^ 9 + cos"^ 6 + sin^ 6 cos^ 9 cos (C — Svr/d) , (67) 

8 V 2 

Pe{xn, Xc) OC (^1 (C “ 3 ^/ 4 )) siu^ 9 cos^ 9 , ( 68 ) 

where 

C = mii/^^. (69) 


It is important to stress the difference between the cases of sufficiently extended wave 
packet source and the hxed space-time point source. In the hrst case the probability of 
detecting neutrino is given by Eq. (|33|) with the phase given by Eq. (^), in the second case 
the probability is given by Eq. (|B^. The origin of this difference lies in the fact that in 
the former case the amplitude (|l|) involves integration over Xc, whereas in the latter case the 
point Xc is hxed. In the asymptotic limit in which Eq. ( [55|) is valid neutrino propagators have 
strong pre-exponential mass dependence that results in the peculiar mass-dependence of the 
probability (|6^) . With propagators in the asymptotic limit (|59|) it can be explicitly demon¬ 
strated that integration of the amplitude (|^ over Xc in the case of monochromatic source 
produces neutrino-mass-dependent factors that cancel out such pre-exponential neutrino- 
mass-dependence of the probability amplitude and also modify the phase of the probability 
amplitude, leading to Eq. (|30|) . In view of the derivation of Eq. (^OD presented in Sec. ^ 
such a demonstration in the main text would be redundant. Therefore, in order to make this 
point clear, we perform it in the Appendix. 

In reality, creation of neutrino cannot occur at a hxed space-time point. However, possible 
creation region might happen to be sufficiently localised by the nature of the source or by the 
experimenter, so that the phase diherences between components of (|6T|) will be well hxed. 
Thus the equations of this section will apply to the situation when neutrino creation region 
is localised in space and time in such a way that 



where by 6[f{xc)] we signify characteristic variation of f{xc) due to variation of Xc over the 

creation region. As the probability oscillation phases mjk\J[xn — Xc)^ are symmetric with 
respect to Xn ^ Xc interchange, this means that the creation region is to be restricted in 
space and time by the oscillation time and length scales in the vicinity of the point Xn- Since 
for small variations 


5 


Tfljk 




^nc 


■ Sxc 



(71) 


these oscillation scales will be determined by the four-momentum diherences pjk = 


xdi with hxed four-vector Xr, 


'^jk ^nc 

detection event x. 


xd which is the neutrino four-velocity u at the 


14 










Within the model of two neutrinos considered above we will have 


Am 

Pi — P 2 = uAm = Pi -, 

mi 


( 72 ) 


where pi is the four-momentum of the heavier neutrino mass specie at the detection point. 
In the case of close neutrino masses, Am -C mi, the oscillation length and time will be given 
by 


— -*■ osr — 


mi 


EryAm 


mi 


(73) 


where L is given by Eq. (|37D and En is the neutrino energy at the detection point. Then the 
equations (|6^, (|66D will be valid as long as 




(yt < 


ml 

K 


(74) 


In the case of mass hierarchy, mi S> m 2 we have 

P 1 -P 2 -P 1 , (75) 

so that 

Lose - Lose ^ , (76) 

where Ei is the energy of the heavier neutrino at the detection point. The equations for the 
probabilities in this case are valid for 


a^, at<E^K (77) 

Due to incoherent distribution of the sources in realistic situations (in a supernova or 
in the Sun) the probabilities (|65| ) - (|68|) will be averaged over the sources and the oscillatory 
terms are likely to disappear. Whenever this is the case, the probabilities will describe the 
so-called global effects of appearance-disappearance of neutrino flavour species, and will look 
like 


Xc) oc ml cos^ 9 + ml sin^ 9 , 
Pe{xn, Xc) oc (ml + mfj sin^ 9 cos^ 9 , 


in the case mj 



S> 1, and 

Pfiixn, Xc) OC sin^ 9 cos"^ 9 , 

sin^ 9 cos^ 9 


Pe{Xn,Xc) OC 


TT 

1 + W 


sin^ 6^ cos^ 0 , 

8 


(78) 

(79) 


(80) 

(81) 


in the case m 2 ij 1 <C mujx^^^, where ( is given by (l6^) . 

To see whether the assumption of a well-localised neutrino source is realistic, let us make 
some estimates for the cases of solar and supernova neutrinos. We take all the data from 
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the book p. In all these cases one has ax -C at, and it is the value of at which will be 
important .[| In the case of solar neutrinos 

at ~ 10“^ cm, En ^ 10 MeV. (82) 


With Ai2 = ml 


mo ~ 


10 ^ eV^ in the case of close neutrino masses the condition 


under which the formulae of this section will apply, will read 


mi > 3eV. (83) 

Since E~^ ~ 10“^^ cm, in the case of mass hierarchy the formulae of this section will not 
work. 

For supernova neutrinos from the core 

at ~ 10“^^ cm, En ^ 100 MeV, (84) 

so that E~^ ~ 10“^^ cm, the condition (|77D will be on the edge of fulhllment and the formulae 
of this section might be applicable. 

For supernova neutrinos from the neutrino sphere 

at ~ 10“® cm, ~ 10 MeV, (85) 


in the case of mass hierarchy the expressions of this section will not be applicable. In the 
case of close neutrino masses for A 12 = ml — m\ ~ 10““^ eV^ we will have the condition 


mi > 0.3 eV, 


( 86 ) 


for which the relevant expressions of this section will apply. 

In connection with the above examples we must note that the probabilities (|6^ ) -(|68| ) will 
refer to neutrinos as they appear from the source. Subsequent neutrino scattering off the 
particles of the solar or supernova media will result in the well-known Mikheyev-Smirnov- 
Wolfenstein effect (see P, 0]), which is not considered in this paper. 


6 Summary 


In this paper we treated the problem of neutrino flavour oscillations by consistently using 
space-time description of the relevant processes of particle creation and subsequent detection. 
We described the EPR-like experiments of detecting neutrino together with the accompany¬ 
ing charged lepton, as well as the standard textbook examples of neutrino flavour space-time 
oscillations, without invoking a priori the notion of particle trajectories. From our analysis 
it is also clear why in fact it is possible to use such a notion. The effective integration region 

^Note that in Q as well as in some other literature tr^ usually stands for the emitted neutrino wave 
packet spread. In this paper both Ux and cr* denote the spread of the neutrino source. Also note that we 
put the speed of light as well as the Planck constant to unity and measure at in units of length. 
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(the region of constructive interference) in the probability amplitude (|^ over the space-time 
point Xc of particle creation is localised around the place determined by particle classical tra¬ 
jectories, and the contribution to the phase of the amplitude comes mainly from the action 
along these trajectories. In the case of wave packet neutrino source our treatment enabled us 
to obtain in a rather simple way the coherence lengths of the oscillations. We also considered 
the case of neutrino source strongly localised in space and time and in this case found out 
dependence of the probability oscillation amplitude on the neutrino masses. 
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Appendix: Alternative derivation of the probability am¬ 
plitude. 


In this Appendix we shall derive the expression for the probability amplitude (|30D in the 
case of monochromatic source, in the limit mn\Jx‘^^ 3> 1, S> 1, using the asymptotic 

expression (^) for the propagators in the coordinate representation. The integral over Xc 
in (|^ will be evaluated in the stationary phase approximation, according to the assumption 
that the main contribution comes from the region of stationary phase of the integrand. This 
phase stems from the propagators and from the source current in (|^ and is given by the 
expression 




^{xi,Xn,Xc) = - 

and it’s stationary point is determined from the condition 


p- Xc 


dx° 


= mixfc + mnx'^c 


-P =0, 


(Al) 


(A2) 


where the notation is used x = x/yx"^. Since xic and Xnc are just four-velocities, respectively, 
of the charged lepton and the neutrino at their respective detection points, the equation ([A^ ) 
expresses the total energy-momentum conservation. Note that for different neutrino masses 
mn = mj the value of Xc determined by Eq. (^) will be different. Let Xj be the solution 
for Xc of Eq. ([A^) with nic = rrij. For further convenience we make notation xi — Xj = Xij, 
Xn — Xj = Xnji xij = Ui, Xnj = Un- The phase d) can be developed in powers oi x = Xc — Xj 
around the stationary point a: = 0 with the result 


^{xi,Xn,Xc) = -Cafsx'^x'^ + 


(A3) 
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where 


(A4) 

(AS) 


Ccif^ Cliui)Q^{ui) Cjii^Un)ai^Un) (q “|“ Cjj) Qaf3 ; 



Dropping the higher order terms in (A3), denoted by dots, we will be interested in the value 
of a Gaussian integral 

(A6) 


d'^xexp (^Cajjx'^x^^ 


This, up to a constant factor, is given by |det{GQ,^}|~^'^^. In terms of the velocities vi and 
respectively, of the charged lepton and the neutrino in the source rest frame the determinant 
is given by 


det{G„/3} = ClCn {ci + Cnf {ui ■ Unf “ 1 = QC^ (q + Cn) 


{Vl + Vn 


(1 - v^) (1 - n2) 


(A7) 


The values of q and c„ given by ( [A5| ) can be easily seen to coincide with the extremal values, 
respectively, of A; and A„ determined by the system of equations (|[), (|5), and given by (|^) . 

Combining all the factors in (|^) together we obtain the expression for the amplitude in 
the limit of —>■ 0 for the prefactors and up to an irrelevant constant as 


DCn D ) (7 ^n) ^ oir\ /l Tt\ TT 

^ I ^ ^ - OolaUnOs [I - ) U^jUja 6 




(ci + Cn) 


Vl + Vr, 


md 


OolaPnOs {m - ^ UljUja , 


(AS) 


where is the value of the phase at the stationary point of Xc = Xj that corresponds to 


neutrino mass specie j. In view of the expression (^) for the phase we see that this value 
is identical to that of Eq. P3|). Therefore the expression ( |A8|) for the amplitude coincides 
with that of Eq. (|30|) . 

Note that the pre-exponential factors in the hnal expression (|A8|) for the amplitude 
remain hnite in the limit of rUn = 0 in spite of the fact that neutrino propagators (^U] ) have 
pre-exponential factors The reason is that the strong neutrino-mass-dependence of 

these factors has been counterbalanced by the neutrino-mass-dependence of the value x\ 


nj 


as well as of the factor |det{Ga/ 3 }|~^'^^, with the determinant given by ( |A7| ). In the case of 
hxed neutrino creation point there is no integration over Xc and factors remain in the 
probability amplitude, Eq. (pT|). 
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